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ABSTRACT 

It  is  shown  that  the  adjoints  of  a  spline  based  approximation  scheme  for 
delay  equations  do  not  converge  strongly. 
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1.  INTRODUCTION 


I 

I 

The  primary  goal  of  this  paper  is  to  illustrate,  by  a  simple  problem,  the 
necessity  of  conducting  a  careful  analysis  of  numerical  schemes  (that  were  de¬ 
veloped  for  open-loop  simulation)  when  these  schemes  are  to  be  applied  to  an 
'  optimization  based  control  design  problem.  It  is  our  feeling  that  many  dis¬ 
tributed  parameter  control  systems  (viscoelastic  structures,  fluid  flow  con-  f 

trol,  etc.)  are  such  that  "standard  finite  element/finite  difference"  schemes 

might  lead  to  numerical  difficulties  in  certain  control  design  methods,  if 
-  < 
care  is  not  taken  to  ensure  that  these  approximation  schemes  have  convergence  j 

properties  essential  for  the  intended  application. 

We  shall  concentrate  on  a  non-convergence  result  for  an  optimal  control 

problem  governed  by  a  delay  differential  equation.  Although  we  feel  that  the 

technical  details  required  to  prove  non-convergence  are  interesting,  we  hope 

that  this  proof  is  not  viewed  as  the  major  contribution  of  the  paper.  Indeed, 

we  hope  that  the  reader  is  motivated  to  think  about  similar  problems  for  more 

complex  distributed  parameter  systems. 

During  the  past  fifteen  years  considerable  attention  has  been  devoted  to  5 

i 

the  construction  of  finite  dimensional  approximations  of  distributed  parameter 
systems.  Much  of  this  work  is  based  on  algorithms  first  developed  primarily 
for  simulation.  However,  it  is  not  clear  that  finite  dimensional  models 
developed  for  (open-loop)  simulations  will  also  be  suitable  for  certain  opti¬ 
mization  based  control  design  techniques.  Moreover,  there  may  be  several 
reasons  that  a  numerical  scheme  developed  specifically  for  simulation  does  not 
perform  well  when  applied  to  a  control  design  problem. 

In  this  note,  we  concentrate  on  the  use  of  the  so-called  spline  scheme 
developed  by  Banks  and  Kappel  [2]  as  an  approximation  method  for  calculating 


optimal  feedback  gains  for  control  systems  governed  by  functional  differential 
equations.  This  problem  is  simple  enough  to  be  addressed  in  a  short  paper  and 
yet  it  still  illustrates  how  a  specific  numerical  scheme  when  combined  with  a 
particular  control  design  approach  can  lead  to  numerical  difficulties. 

In  order  to  motivate  the  paper  and  to  describe  the  main  technical  con¬ 
tribution  of  the  paper,  we  first  review  some  known  results. 

Let  H  and  U  be  real  Hilbert  spaces  and  S(t):  H  H  denote  a  Cq- 
semigroup  of  bounded  linear  operators  with  generator  A.  We  assume  that 
B:  U  +  H,  Q:  H  -*■  H  and  R:  U  +  U  are  bounded  linear  operators  with  Q 

and  R  self-adjoint,  non-negative  and  R  satisfies  R  ^  ml  >  0.  The  (LQ) 
optimal  control  problem  is  to  minimize 

T 

(1.1)  J(u)  =  /  [<Qz(s),  z(s)>  +  <Ru(s) ,u(s)>]ds 

0 

where  z(t)  is  defined  by 


t 

(1.2)  z(t)  *  S(t)zn  +  J  S(t-s)Bu(s)ds 

U  0 

for  0  £  t  <  T  and  given  Zq  €  H.  If  T  -  +»,  then  one  has  the  linear 
quadratic  regulator  problem. 

Assume  now  that  there  exists  a  sequence  of  Cg-semigroups  SN(t)  on  H 
and  positive  constants  M,  6  such  that 


(1.3) 


a  SN(t ) I  <  Me6t,  t  >  0,  N  >  1 


and 


-3- 


(1.4) 


SN(t)  ♦  S(t)  strongly  as  N  -*■  +» 


and  the  convergence  in  (1.4)  is  uniform  in  t  on  bounded  intervals.  Denote 

by  the  generator  of  S^(t)  and  also  assume  the  existence  of  bounded 

linear  operators  B^:  U  -*•  H,  QN:  H  ■*  H  with  QN  self-adjoint  and 
N 

Q  >0  satisfying 


(1.5) 


ibn  -  B«  -»■  o,  «qn  -  Q1  +  0. 


Note  the  assumption  of  uniform  convergence  (1.5)  is  stronger  than  required  by 
Gibson  (see  page  113  in  [4]).  It  is  well  known  that  the  optimal  control  (if 
it  exists)  is  given  by  state  feedback 


(1.6) 


u  (t)  -  -R  H(t)z  (t) 


where  the  bounded  self-adjoint  operator  II(»)  satisfies  the  Riccati 
(integral)  equation 


(1.7) 


n(t)z  -  /  S*(s  -  t) [Q  -  n(s)BR_1B*n(s) ]S(s-t)ds. 
t 


Let  H  (t)  be  the  solution  to  the  "approximating"  Riccati  equation 


(1.8)  nN(t)z  *  /  SN*(s-t)[QN  -  nN(s)BNR"1BN*nN(s)]SN(s-t)ds 


and  observe  that  (1.8)  would  determine  the  sub-optimal  gains  if  one  used  the 
approximating  system  (A^,B^)  with  weights  Q1*. 


u 


The  following  theorem  Is  a  direct  consequent  of  Gibson's  (more  general) 
results  (Theorem  6.1  and  Theorem  6.2  in  [4]). 

Theorem  1.1:  If  conditions  (1.1)  -  (1.5)  hold,  then  for  0  £  t  <_  T, 

N 

II  (t)  converges  weakly  to  n(t)  and  the  convergence  is  uniform  on 
[0 ,T] .  If  in  addition 

(1.9)  SN*(t)  -»■  S*(t)  strongly, 

then  nN(t)  converges  strongly  to  n(t)  and  the  convergence  is  uniform 
for  t  €  [0,T] . 

Corollary  1.2:  Let^  K(t)  «  R-1B*n(t)  and  KN(t)  =  R_1BN*nN(t) 

denote  the  feedback  gain  operators.  Assume  that  U  *  IP  (i.e.,  is  finite 
dimensional).  If  nN(t)  converges  strongly  to  n(t),  then  as  N  -*■  +», 

(1.10)  IKN(t)  -  K(t)l  -*•  0. 

The  main  point  to  be  emphasized  is  that  if  the  control  space  is  finite 
dimensional,  then  uniform  convergence  is  assured  provided  the  numerical  scheme 
is  stable  and  consistent  and  (1.9)  holds.  If  one  is  concerned  only  with  simu¬ 
lation,  then  stability  and  consistency  is  sufficient  for  most  numerical  prob¬ 
lems.  Moreover,  it  can  be  shown  that  many  standard  numerical  schemes  de¬ 
veloped  for  simulation  of  self-adjoint  partial  differential  equations  do 
satisfy  (1.9).  Therefore,  it  is  not  surprising  that  until  Gibson  "needed" 


(1.9)  to  establish  the  uniform  convergence  of  optimal  gain  operators,  the 


question  of  whether  a  numerical  scheme  satisfied  (1.9)  received  little  atten¬ 
tion.  Indeed,  even  after  Gibson  published  his  result  it  was  still  not  obvious 
that  condition  (1.9)  was  anything  more  than  a  technical  assumption  needed  in 
Gibson's  proof. 

In  [3]  Banks,  Ito,  and  Rosen  applied  a  convergent  spline  based  scheme  to 
an  optimal  control  problem  governed  by  a  delay  differential  equation.  The 
numerical  results  in  [3]  seemed  to  show  that  KN(t)  did  not  converge  uniform¬ 
ly  to  K(t)  and  these  numerical  results  have  often  been  used  as  evidence  that 
(1.9)  did  not  hold  for  this  particular  scheme.  Moreover,  several  new  schemes 
have  since  been  generated  specifically  to  ensure  (1.3)  -  (1.4)  and  (1.9)  are 
valid.  Still,  it  was  not  known  if  condition  (1.9)  held  for  the  spline  scheme 
used  in  [3].  We  shall  provide  a  proof  that  (1.9)  fails  for  this  scheme.  We 
also  show  that  this  spline  scheme  is  stable  and  consistent  to  A  on  a  dense 
subset  of  D(A*) . 

2.  SPLINE  APPROXIMATIONS  OF  HEREDITARY  SYSTEMS 

Consider  the  delay  differential  equation 

0 

(2.1)  x(t)  *  AQx(t)  +  AjX(t  -  r)  +  /  A(s)x(t  +  s)ds 

-r 

with  initial  data 

(2.2)  x(0)  ■  r\ ;  x(s)  -  $(s),  -r  £  s  <  0, 

where  x(t)  I?  and  the  elements  of  A(«)  are  square  integrable  on 
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[-r,0].  It  is  well  known  that  (e.g.,  [1])  for  n  €  HP  and 

4>  €  L  (— r , 0 ;  if1)  there  exists  a  unique  solution  of  (2.1)  -  (2.2) 

x:  [-r,+°0  lP  such  that  x£W^’^(0,T;  I?)  for  all  T  >  0.  If  one 

defines  the  solution  map  S(t),  t  >  0  on  the  product  space 

Z  =  I?1  x  L2(-r,0;  l?1)  by 


(2.3) 


S(t)(n  ,♦(•))  =  (x(t),  x(t+.)) 


where  x  is  the  solution  to  (2.1)  -  (2.2),  then  {S(t)}^Q  is  a  strongly 
continuous  semigroup  (i.e.,  C0-semigroup)  on  Z.  The  infinitesimal 
generator  A  is  the  operator  defined  on  the  domain 


(2.4) 


D(A)  =  {0l  ,$(•))€  Z |<j>(»)  €  W1,2(-r,0;  IP),  $(0)  =  n} 


(2.5) 


A(n  ,<t> ( •  ) )  *  (Aqti  +  (-r)  +  /  A(s)<(i(s)ds,  $(•)). 


The  adjoint  operator  A*  generates  the  adjoint  semigroup  S*(t)  and  it  is 
easy  to  show  that  (see  [4,7]) 


(2.6) 


D(A*)  =  {(5  ,<!>)€  Z|ij>e  W1»2(-r,0;  If1),  ^(-r)  =  A^} 


and  for  (£  ,iji)  €  D(A  ) 


(2.7) 


A" U  ,*)  -  (1<(0)  +  aJs,  [AT(.)5  -  $(•)]). 
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As  In  [2],  we  define  the  linear  spline  based  approximation  for  S(t). 


Let 

iRNlN 

{Bi}i=0 

denote 

the  usual  linear 

B-splines  defined  on  the  interval 

f-r,0] 

by 

N  .  N  . 

17  <s  -  W' 

,  r  N  N. 

se[Ti+l>Ti] 

(2.8) 

B^(s)  =  < 

,  r  N  N  , 
s  e  [Ti,  T±_1  ] 

1 

otherwise. 

where 


N 


Ti  =  -ir/N,  i  =  0,1, tn+1  =  -r  and  t_1  =  0 


For  each 


N  =  1 , 2 ,  • • •  let  ZN  denote  the  linear  subspace  of  Z  defined  by 


(2.9) 


ZN  =  {z€Z|z  =  l  a^CB^CO)  ,B^(-  )),  a^eS?} 
k=0 


and  let  PN  denote  the  orthogonal  projection  of  Z  onto  ZN.  This  subspace 
can  be  identified  with  jgi(N+l)  ^  the  prolongation  iN:  ip(N+l)  +  z 

defined  by 


N 


(2.10) 


A  ■  <a0-  I  WN<-» 

k=0 


where  a  =  (aj,  aj,*  •  •  ,a^)^€  .  The  space  |^(N+0  is  normed  with  the 

induced  inner  product 


(2.11) 


<a,b>N  =  aTQNb, 


where  a,  b  €  lfl(N+1)  and  QN  is  defined  by 


i*  ••  '■*  w 


■BMunnasK 


'  VI 

*  VI 


i 

i 


•u*’ 

$*; 

’,'s 


m 


*,y 

vl.i 

ift 

iiS 

y 


i 


?s 

$ 

$ 

!®| 

jk 

'<} 

M* 

8 

M 


■5 


•J  <| 

*<l 

a 

.■I 

*i? 


:: 


(2.12)  QN  =  0 
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The  adjoint  operator  [iN]*:  Z  ■*  #(N+1)  is  given  by 


(2.13) 


n  +  ♦  , 


[iN]*(n  ,4> (• ) )  «  [QNf 1  ?1 


N  /•  N 

where  <t> ^  =  J  4>(s)B^(s)ds.  Moreover,  it  is  easy  to  show  that 
-r 

{  [iN]  iN  =  I  the  identity  on 


(2.14) 


iN[iV  =  PN 


and  for  z,  w€  Zr 


(2.15) 


<z,w>z  =  <[iN]*z,  [iN]*w  >N. 


In  order  to  construct  the  standard  Galerkin  approximation  of  A,  we  note 
that  ZN  C_  D(A)  and  define  AN  by  A^  =  P^AP^.  Observe  that  AN  (and 
hence  [An]*)  is  continuous  and  although  P%  c_  D(A),  PN  does  not  map  all 
of  Z  into  D(A*).  It  is  shown  in  [2]  that 


an  -  iVrVuV 


(2.16) 


where 


(2.17) 


.N  aN  4N 

-1  .  o" 

A0  A1  — 

I  0 

.  1 

•  • 

0 

+  2 

•  *  0  *-I 

0  *1-1 

and 

0 

Aq  =  Aq  +  /  A(s)Bg(s)ds 

-r 

N  0  N 

\m  f  A(s)B^(s)ds,  1  <  k  <  N-l 
-r 

A^  =  A1  +  /  A(s)B^(s)ds. 


If  SN(t)  denotes  the  Cg-semigroup  generated  by  AN,  then  it  is  shown  in  [2] 
that  for  each  z€Z, 


(2.18)  ISN(t)z  -  S(t)zK  0, 

where  the  convergence  is  uniform  in  t  on  bounded  intervals.  The  convergence 

(2.18)  was  established  by  proving  that  the  Galerkin  approximations  AN 
satisfy  the  Trotter-Kato  Theorem  [6]  and  hence  provided  a  stable  and  consis¬ 
tent  approximation  scheme  for  A.  We  shall  prove  that  the  above  convergence 
statement  does  not  hold  for  the  sequence  of  adjoint  semigroups  [SNj*(t). 
Moreover,  this  convergence  fails  even  though  the  adjoint  operators  [AN]*  are 
stable  and  consistent  to  A*  on  a  dense  subset  of  Z  (i.e.,  there  is  a  dense 
set  C  C  D(A*)  c  Z  such  that  [A^]*z  +  A*z  for  all  z €  C) . 


3.  CONVERGENCE  OF  THE  ADJOINT  GENERATORS 


We  shall  follow  the  approach  given  in  [5]  and  consider  the  n(N+l)  x  n 
matrix 

N  r^N.-l,.  _  tNT 

Y+  =  (Q  ]  (0 . 0,1)  . 

Define  the  operators  6  +  :  IFf1  by 

6+x  =  4>x(*  )) 

where 

♦x(0  =  [Bq(» )i,  bJ(Oi,  bJ(.)i]yJx. 

It  follows  that  for  all  (n  ,4> )  =  z€Z  and  xflP 

rf.N  *-N  .T.Nf.N,*  N  r  N  .T-N.-N,* 
l l i  1  6+x]  Q  [i  ]  P  z  =  [y+x]  Q  [i  ]  z 

T  N  N  N  * 

=  xT[B"(-r)I,...,Bj(-r)I][iN]  z, 

and  (2.14)  -  (2.15)  implies 

(3.1)  <6^x,z>z  =  xT[$N(-r)], 

where  PNz  =  P**(n,$)  =  (<J> ^ (0 )  ,$^(«  ) )  e  Z^.  Furthermore,  if  ^^in  denotes 
the  smallest  eigenvalue  of  Q  ,  then  (r/6N)  £  X  ^ 


and  it  follows  from 


We  also  need  the  following  representation 


N  *  N 

Lena  3.1:  The  operators  [A3  :  Z  Z  are  given  by 

[AN]  =  pV<0)  +  aJ/(0),  [AT(.HN(0)  -  4iN(»)]) 

r N /  ,T  N,.v  N ,  N« 

+  6  +  (A1i(;  (0)  -  ip  (-r)), 


where  PN(5  ,i|0  =  C^N(0)  ,4>N(*  )  ) . 

Proof :  Assume  that  z  =  (r),<f>)  and  w  =  (£,<|>)  belong  to  Z  and 
let  P^z  =  (4.N(0)  ,4>N(* ) )  and  P^w  =  (ij)^(O)  ,i|>^(*))  denote  the  orthogonal 
projections  of  z  and  w,  respectively.  The  identity  (3.1)  implies  that 

<0|>N(0)  +  aJAo),  (at(.)*n(0)  -  ♦NC-)1).  pNz>z 

+  <6^(A^N(0)  -  ^N(-r)),z>z 

=  [/(0)]T[A0)]  +  [/(0)]TA0UN(0))  -  /  <^N(s)  ,4>  N(  s )  >ds 

-r 

+  [Ao)]T  f°  A(sHN(s)ds  +  [^N(0)]TA1[«N(-r)] 

-r 

.  N,  v,Tr  N, 

-  [i|»  (-r)l  [<fr  C-r)  1  - 

Integrating  by  parts,  the  boundary  terms  cancel  with  the  first  and  last  term 
in  the  sum.  Therefore,  it  follows  that  if  [A^]  is  defined  as  above,  then 
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,  N  *  N  N  N 

<[A  ]  w,z>  =  <P  w,AP  z>  =  <w,A  z> 


and  this  completes  the  proof. 

It  should  be  noted  that  [A^]*w  =  P^A  P^w  if  and  only  if 
Ay(0)  -  xpN (— r )  =  0,  i.e.,  if  and  only  if  PNwCD(A*).  Also,  if 


D  =  D(A)H  D(A*) 


and 


I 


C  =  {(£  ,i|i)£  D|*  e  C2(-r ,0;  I?1)} 


then  D  and  C  are  dense  in  Z.  Moreover,  we  have  the  following  convergence 
result. 


Lemaa  3.2:  C  and  D  are  defined  as  above,  then 


(a)  [A^]*w  ♦  A*w  for  all  w€C 

and 

(b)  for  all  X  €  ®,  (XI  -  A*)D  is  not  dense  in  Z. 

Proof :  Let  w €  C  C  D(A)  O  D(A*) .  Note  that  w  *  0|>(0)  ,i|>(» ))  where 
(— r )  =  A^ip ( 0 )  and  PNw  -  (^N(0)  ,^N(*  ) )  €  D(A)  .  It  follows  from  (2.7)  and 
Lemma  3.1  that 


N.  *  *  m  *  * 

[A  ]  w  -  A  wl  <  IP  A  w  -  A  wH 


+  .(*%)  +  Aq^P  N  ( 0 ) ,  [AT(*  )ipN(0)  -  $N(*)])  -  (*(0)  +  A^(0)  ,  [Aq(»  )i(»(0)  -  ^(»)]) 


Since  w=  (i|;(0)  »  <K*  ))  €  C  C  D(A)  Pi  D(A*)  and  £  C2(-r,0;  1?),  it 

follows  from  standard  estimates  on  interpolating  splines  (see  equations  (4.1) 
-  (4.3)  in  [2])  that 

|/(0)  -  *(0)|  <  0(1/N2) 

|/(s)  -  «(s)|  <  0(1/N),  -r  <  s  <  0 
and 

»/(•)  -  ,J(.)j  <_  0(  1/N). 

The  first  term  +  0  since  w€D(A*)  and  pNz  -*■  z  for  all  z€Z.  The 

second  term  is  estimated  by 

f”  <  |/(0)  +  A;/(0)  -  <J»(0)  -  A>(0)|  +  BAT(.)(/(0)  -  *(0))l 


+  */(•  )  -  ^(»  )# 


<  (1  +  | Aq|  +  IA(-)l)|/(0)  -  iK0)|  +  0(1/N) . 


Therefore, 


„N 


0 


as 


N  » 


Applying  (3.2),  the  last  term  is  estimated 
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<  (6N/r)1/2(|A1|  |*N(0)  -  ^ CO)  |  +  |*(-r)  -*Vr)|), 

N 

and  hence  +  0  which  establishes  part  (a). 

Turning  to  part  (b),  let  w  *  (i|r(0)  ,i|i)  £  D,  5 ,  x  £  if1  and  X  €  *.  A 
straightforward  calculation  yields  that 

<(5»e*  x),  [XI  -  A  ](<|>(0),i|i(«  ))> 

=  <|iT(0)[(X  -  1)1  -  Aq]?  +  4iT(0)[I  -  e“XrAj  -  /°  eXsAT(s)ds]x. 

-r 

If  (X  -  l)€o(A0)  and  eXr£  o(Aj  +/  e^X+r  ^8AT(s)ds) ,  then  there 

-r 

exist  5*0  and  x  *  0  such  that 


<(5  ,eX*x)  ,  [XI  -  A* } (4» (0 )  ))>  -  0. 


If  (X  -  1)  4.  cx(Aq),  let  x  *  0  and  define  5  by 


5  -  -[(X  -  1)1  -  A0]"l(I  -  e_XrA1  -  /  eXsAT(s)ds)x, 

-r 


or  if 


eXr  t  o(Aj  +  j  e^X+r^SAT(s)ds) ,  let  5*0  and  define  x  by 
-r 


x  =  -[I  -  e"XrAj  -  /  eXsAT(s)ds]_l [ (X  -  1)1  -  AQ]5. 

-r 


In  any  case,  there  always  exists  an  element  z  *  (5  ,<!»(•))  ■  (5,eX*x)  *  0, 
such  that 

* 

<z,  (XI  -  A  )w>  *  0 
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for  all  w  *  (ij>(0),  4>(*))eD.  Consequently,  (XI  -  A*)D  is  not  dense  in 
Z. 

Since  «[SN]*(t)l  =  ISN(t)l  _<  eut  for  some  uj  independent  on  N 
(see  (2)),  the  existence  of  a  set  E  C  D(A*)  and  a  X  €  F  such  that 
(XI  -  A*)E  is  dense  in  Z  and  [A^]*w  *  A*w  for  w  €  E  would  imply 

strong  convergence  of  the  semigroups  [SnJ  (t)  to  S  (t)  ([6],  III,  Th. 

4.5).  Although  [AN]*w  +  A*w  on  the  dense  set  C,  (XI  -  A*)C  is  not 
dense  in  Z.  We  shall  establish  that  there  does  not  exist  a  set  E  c;  D(A*) 
with  the  above  properties.  In  fact,  we  shall  show  that  [S^]*(t)  does  not 
converge  strongly  to  S*(t). 

4.  CALCULATION  OF  !*( [AHJ*)“1PH 

In  this  section,  we  present  several  technical  lemmas  that  will  be  needed 
in  Section  5.  The  proof  of  Lemma  4.1  is  straightforward  and  hence  omitted. 
Lemma  4.3  can  be  found  in  [4,7]. 

Lean  4.1s  The  operator  [AN]*  can  be  decomposed  as 

(4.D  uV  -  iV^nmiN]*. 


Lemma  4.2:  Assume  that  is  invertible.  If  z  »  (n  ,$(•))€  Z, 

then  PN( [A^]*)”^P^z  exists  and 


(4.2) 


pN([AN]*)-lpNz  .  iNaN( 
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where  is  the  solution  of 


(4.3) 


[HN]TaN  = 


n  + 

,  N 


bN 

N 


N 


Proof :  Consider  the  equation  A^*  =  PNz  for  w  €  Z^.  By  Lemma  4.1  and 


from  (2.14) 


. N/nNv -1UNT  . N*  ,N,N* 
i(Q)  H  i  w  =>  i  i  z. 


Multiplying  with  i^  and  from  (2.13)  and  (2.14) 


UNT  .N* 

H  i  w 


n  +  $ 
,  N 


an 


N 


M  M  M*  M  M 

The  lemma  now  follows  from  the  fact  that  w  ®  P  w  =  I  i  w  =*  i  a  . 


4.3:  If  A  =  A 


Q  +  Aj  +  /  A(s)ds,  then  Oep(A  )  if  and 


-r 


only  if 

(4.4) 


A  is  Invertible  and  (5  *♦(•))  “  (A  )  (n,<K*))€D(A  )  is  given 


C  *  [A1]-1 (n  +  /  ♦(s)ds) 


-r 


S  8 

(4.5)  <Hs)  -  (aJ  +  /  AT(0)d9)5  -  /  41  (0  )d0 

-r  -r 


Consider  equation  (4.3).  It  follows  from  (2.17)  that  (4.3)  is  equivalent 


to  the  system 


! 


•S' 

;S 

!«» 


*«ir‘ 

id 


V, 

,y: 

'Si 


t!*Sj 

$ 

ys 


•.* 

■> 

:■> 


I 


(4.6) 


(4.7) 


(4.8) 


1  /-  N  ,N  \  j.  r*NJ  N  ,N 

2  (  k+1  "  ak-l}  +  lAk]  a0  “  *k  1  <  k  <  N_1 


1  / _N  ,  N  w  f .NiT  N  _ 

”  1  (aN  +  aN-l)  +  [V  a0  ' 


1  ,  N  .  Nw  .  .N,T  N  .  XN 

2  (a0  +  31)  +  {A0]  a0  *  n  +  ^0* 


01  0 


A  straightforward  induction  argument  yields 


(4.9) 


N  _  N  o  V  /r*N  ,T  N  N  .  .  „ 

a2k  a0  ~  2  Z.  ^  A2i— 1  a0  ^2i-P  *  ^  2^  —  ^ 


(4.10) 


a2k+l  *  a?  -  2  l  (fA2i]Ta0  "  +2i)*  lf  2k  +  1  <  N* 
1=1 


Thus,  it  follows  from  (4.7)  that 


(4.11) 


1  .  »  N.  .  ?  ,.N.T  N  ?  ,» 

-J  <a0  '  *1>  +  I  a0  '  l,  V 

k=l  k*l 


Moreover,  (4.8)  and  (4.11)  imply  that 


where 


N  NTH  N  0 

I  a  -  r,  +  l  -  h  +  /  (Ks)ds 

c*0  j-0  J  -r 


N  0 

I  f )T  -  [A.  +  A.  +  /  A(s)ds]T  -  AT. 

k*0  K  1  -r 


*  T 

If  one  assumes  that  0ep(A  ),  then  by  Lemma  4.3  it  follows  that  A  is 
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invertible  and 


(4.12) 


N  —  T  ® 

aQ  =  A  (n  +  /  4>(s)ds), 

-r 


N 


Observe  that  (4.12)  implies  that  a^  is  independent  of  N  and  by  (4.4)  it 
follows  that  ao  =  ^  where  (5,4<(0)  *  [A*]-1  (n  (•  ) )  .  Equation  (4.8) 


yields  the  identity 


(4.13) 


N  N  .  .  ?  /r*N,T  N  Nx 

31  “  “ao  +  2  (IV  a0  '  *k} 

k=l 


and  hence  it  now  follows  from  (4.10)  that 


N 

*2k+l 


(4.14) 


N  .  -  r  Man,t  N  Nx  ,  r  m.N  ,TN  N 

0  +  2  ki1  (lAk1  a0  *k>  2  ([A2i]  a0  *2i> 

N  .  9  v  .N.  .  r  ,,.N  ,T  N  N 

_a0  +  2  l  ( [ A  ]  a  —  ♦  ^  +  2  2.  '  ^A2i-1  a0  ~  *21-1^' 

j-2k+l  J  u  3  i-1  1  u 


Note  that  for  1  <  k  <  N  -  1 


kl!|2  "  1/  ♦(s)Bk(s)ds|2  <  (^J)  /  U(s)  1 2ds 


-r 


N 

rk+l 


and  similarly 


knl2  1  (3n}  /  !♦  (»)  |  2ds 

N 

T1 


tn 

1 2  <  (^n-)  /  !♦  (») l2(l8' 

-r 


Therefore,  (4.9)  and  (4.10)  imply 


•a2kl  i  I aQ I  ^  +  2lAil  >  +  2 /¥7T  (IAI  |aj|  +  1*1  ) 

and 

I a2k+i I  <  !a?l  +  2IAJ  lao!  +  2/771  (,AtL2faol  + 
respectively.  The  identity  (4.12)  yields  the  estimate 

laol  i  lA~Tl  (1  +  *^)Kn,*(*))iz 

and  (4.13)  leads  to  the  bound 

|a“|  -  |-(1  +  2[A^]T)aJ  +  2(n  +  ♦  ") | 

<  |aj|  +  2 |n I  +  2/773N  0 A«L^  |*J|  +  i*il^). 

Combining  these  estimates  one  has  that  for  0  _<  k  _<  N 

|«JJ|  <  Ml  (n  ,♦  (• ) )■  z 

where  M  ^  0  is  independent  of  N.  An  application  of  Lemma  4.2  yields  the 
estimate 

IPN[AN*]_1PNZIz  -  liNaN«  <  (/T+r)M  lzlz 
for  all  z  €  Z.  We  summarize  these  results  in  the  following  theorem. 
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Theoren  4.4:  If  0€p(A),  then 


0Nr.N*1-l0N(,  A/  ..  N  N 
P  [A  ]  P  (n,4>(*))  =  i  a 


where  aN  =  col(  [a^ ]T,*  •  •  ,  [a^]T)  is  given  by 


a!?  =  A  T(n  +  /  <fr(s)ds). 


-r 


and  for  0  <  2k  <  N 


N  N  .  ?  , ,  .N  ,T  N  N  . 
3 2k  a0  “  2  E  A2i-1  30  "  ^2i-I ^ 

i=l 


while  for  1  <  2k+L  <  N 


N  N  .  .  ?  /r*NiT  N  .  >>  r  /  r  ,T  N  N 

*2k+l  =  'a0  +  2  Jk+1  (  Aj  a0  ~  ♦  j )  +  2  ([A^.jl  aQ  -  *2l_ 


M  M*  _  1  M 

Moreover,  IP  [A  ]  P  is  uniformly  bounded  in  N, 


Let  Pave  denote  the  "averaging"  orthogonal  projection  on  Z 
by 

tn 

PAVE I  "  </  ♦(•>d.)X  „  N  ), 

i-1  tn 

where  Xj  denotes  the  characteristic  function  for  the  interval  I 


defined 


(see  [1] 


for  details) 
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Corollary  4.5:  There  exist  a  costant  K  £  0  (Independent  of  N  )_ 
such  that  for  all  z  €  Z 


1DN  ,DNr  .N*,  — 1„N  v  r  a* i-l  „  ✓  ,,  „ 

,PAVE^P  ’  P  z)  -  [A  ]  zll  £  (^.)  Ilzll. 


Proof:  A  direct  application  of  Theorem  4.4  yields  the  identities 


N 


|  (a*..,  +  a".)  =  l  ([A^]TC  -  /),  0  <  2k  <  N-l 


2  '  2k+l  2k' 


i=2k+l 


and 


N 


t  <4-i  *  4>  ■  j2k  -  ♦;>. 2  < *  <  *■ 


*i-i, 


On  the  other  hand,  since  (5  ( •  ) )  =  [A  ]  1  (n  ,<)>(•))  •  Lemma  4.3  implies  that 


5  -  a^  -  A  T(n  +  /  <j>(s)ds) 

-r 


and 


-r 


—  s 

iji(s)  -  (a[  +  /  AT(e)de)e  -  /  <i> (e  )de 

-r 


Therefore,  if  s ^  =  (Tj  +  Tj_j)/2,  then 


l(a2k+l  +  a2k)/2  **s2k+I^ 


N 


2k+l 


I  /  (®2k+l ^ “  l)(AT(s)£  -  <^(s))ds 


N 


2k+l 

N 


2k 


+  /  (B2k+1(8)(A1(8)5  -  $(s))ds| 


N 


2k+l 
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N  N 

T„,  T, 


<  (r//12  N)  {(/  | A(s)  |  2ds)1/2  |C  |  +  (/  ^  |<fr ( s)  |  2ds) l/2}  , 


N 
2k+l 


and  similarly. 


N 

r2k+l 


i<4  +  4-,)/2-*<4>i 


2k- 1 


N 
2k- 1 


<  (r//!2  N)  {(  /  |A(s)|2ds)1/2|c|  +  (  /  |4»(s)  |2ds)1/2> 


2k 


2k 


It  follows  that 


IPN  iNa 

ifave  1  a 


J,  N  N 

j-1  J 


<  (r//12N)  (I Al  |c|  +  l*lT  ). 

l2  l2 


Since  \p  €  W1 ,2(-r,0;  I?)  and 


$(s)  -  Ax(s X  -  *(s)  -  AT(s)A  T(n  +  /  ♦(e)de)  -  *(s), 


-r 


there  exists  a  constant  K  such  that 


- 


l  *<8?>x  N  N  »L  < 

J-1  j  [Tj  2 


(K/N)Kn  ,♦(•))!, 


This  estimate  combined  with  the  previous  inequality  establishes  the  proof. 

We  turn  now  to  providing  a  proof  that  the  approximating  adjoint  semi¬ 
groups  constructed  above  do  not  converge  strongly  to  the  adjoint  semigroup 
generated  by  A*. 


S.  NON-STRONG  CONVERGENCE 

Let  A,  A*.  An,  and  AN*  be  defined  by  (2. A)  -  (2.5),  (2.6)  -  (2.7), 
(2.16)  and  (A.l),  respectively.  The  corresponding  semigroups  will  be  denoted 
by  S(t),  S*(t),  SN(t),  and  SN*(t).  Recall  that  for  z€Z  (see  [2]) 

(5.1)  <ANz,z>  =  <APNz,PNz>  <  o)8PNzS2  <  uj«zB2 
where  w  =  (1  +  2 1  Aq  |  +  |Aj|2+  2D  AS  ^  )/2  and  for  t  0 

(5.2)  ISN(t)«  <  eu>t,  HSN*(t)J  eut:. 

The  following  result  is  a  special  case  of  Theorem  4.2  of  Chapter  3  in  [6], 

Theorem  5.1:  The  following  are  equivalent: 

(a)  For  every  z  €  Z  and  A  €  t  with  Re  A  >  at 

N*  -1  *  -1 

(A I  -  A  )  z  *  (AX  -  A  )  z,  as  N  ♦  «. 

(b)  For  every  z€  Z  and  t  >  0 

SN  (t)z  +  S  (t)z,  as  N  + 

the  convergence  being  uniform  in  t  on  bounded  intervals. 


We  shall  also  need  the  following  technical  lemmas 
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and 


Lemma  5.2:  Suppose  the  condition  (b)  of  Theorem  5.1  holds,  X€p(A* 
II (XI  -  is  uniformly  bounded  in  N.  Then  for  every  z€Z 


PN(AI  -  AN*)-1PNz  +  (A I  -  A*)-1Z. 


Proof:  From  Theorem  5.1,  for  Aq  >  u>  and  z€  Z 


pNv  - 


AN*)_1PNz  +  (XQI  -  A*)_1z, 


as  N  *  <*>. 


Note  that  for  z  e  Z 


PN(XI  -  aN*)~^P^z  -  PN(XqI  -  a^*)-1PNz 


,\Ar  aN*n-I_N,,  T  .N*v  — I  N 

■  Uq  ~  A)P  (XI  -  A  )  P  (XQI  -  A  )  P  z. 


and  similarly 


(XI  -  A*)  !(I  +  (X  -  X0)(A0I  -  A*)"1)  -  (XQI  -  A*)'1, 


’X-l. 


Hence,  if  w  =  z  +  (A  -  Aq)(AqI  -  A  )  z,  then 


PN(XI  -  *N*)-‘pNv  -  (XI  -  aV'w 


(X0  -  X)p"CAI  -  an*)-‘pn((x0i  -  asV‘pns 


<X0I  -  A*)_1z) 


+  p"(X0I  -  a'VVz 


,  *  -1 
(XQI  -  A  )  Az. 


This  implies  that  for  all  z  €  Z 


N,  t  ,N*  -1  N 
P  (XI  -  A  )  P  w 


*  -1 

(XI  -  A  )  w. 


★  “1 

But  since  I  +  (X  -  X Q) (X QI  -  A  )  is  onto,  the  above  statement  holds  for 

all  w  €  Z. 

We  turn  now  to  a  special  case  where  AQ  =  AQ  +  al,  a  £ 

^  ~  ~N* 

A(s)  =  0,  -r  _<  s  _<  0,  A  =  Aq  +  Aj  and  denote  by  A  ,  A  and 

S*(t),  SN*(t)  the  corresponding  infinitesimal  generators  and  semigroups. 

N* 

Lena  5.3:  If  statement  (b)  of  Theorem  5.1  holds  for  S  (t),  then  it 
holds  for  3N*(t). 


Proof:  Note  that 


A 


*  ~M* 

A  +  E  and  A 


N*  N  N 
A  +  P  EP 


where  E:  Z  ■>  Z  is  the  bounded  linear  operator  defined  by 


T 

e(p  ,$(•))  =  (an, -A  (On). 

It  follows  from  (5.2)  that  II  ( t ) 0  <  eut  where 

Consequently,  if  X  >  w ,  then 


u>  =  u>  +  |a  |  +  8  AIT  . 

L2 


J\ 


<*!i 


(5.4)  (XI  -  1  )~L  =  (XI  -  A  )_1  +  (XI  -  A  )_1E(XI  -  A  )  1 , 


Theorem  5.1  implies  that  for  all  z€Z 


(XI  -  AN*)_1z  ■>  (XI  -  A*)_1z 


and  since  the  ranks  of  E  and  E  are  finite  it  follows  that 


N  N 

IIP  EP  -  EK  0. 


It  now  follows  from  (5.3)  -  (5.4)  that 


(XI  -  AN*)_1z  >  (XI  -  A*)  Yz 


for  all  z€Z  and  this  completes  the  proof. 


By  Lemma  5.3,  without  loss  of  generality,  one  can  assume  that  A(»)  =  0 
and  A  =  AQ  +  A^  is  invertible  in  what  follows.  We  will  show  that  there 
exists  an  element  z  C  Z  such  that  P^(A^  )  ^P^z  does  not  converge  to 
(A  )  *z.  First  we  consider  the  case  when  A^  is  not  the  indentity.  From 
Lemma  4.3,  if  (£»♦(•))  =  [A*]  ^(n,0)  where  0  t  p  €  l?1,  then 


5  =  A  and  ijj(s)  =  aJ$  ,  -r  <  s  <  0. 


*  *V*Vm  t 


I*,’ 


•  ,!  *';*,*■-,*»  •  •  *  .*  V  \*  ■»/.»  V  *  \»  V  ,t  ,  1  V  ’  •  '  •  V  '  -  ’  •  ! 
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Applying  Theorem  4.4  we  obtain 


M. 


where 


N, .N*,-lnN,  AS  . N  N 
P  {A  ]  P  (n,0)  =  i  a 


ajk  =  C .  0  <  2k  <  N 


a2k+l  =  +  2A15  »  0  <  2k  +  1  <  N* 


For  illustration,  we  have  the  following  picture  for  the  case  N  -  4  and  n  =  1 


A 

/  \  / 
/  \  / 

/  \  / 


\  /  \  t 

V'.  ■/ - - - v—  i  — 


\  / 
V 


\  / 
V 


T 

Al5 


-C  +  2Aj5 


*  1  N  N 

where  the  solid  line  stands  for  (A  )-1(n,0)  and  the  dashed  line  for  i  a  . 

Since  a!Jk  -  A*5  =  ~  a^R+1  =  ?  -  aJ?  -  e  C  1?  (independent  of  N),  it  is 

easy  to  show  that 


r/2N 


(5.5)  «PN(AN*)_1PN(r,,0)  -  (A*)_1(n,0)ll2  =*  2N  /  |H  s  e  1 2ds  =  e  | 2  *  0. 


Next  we  consider  the  case  At  -  I.  bet  *(s)  -  x  *  0  (constant  vector  in 
Rn).  Then,  from  Lemma  4.3  (£ ,  !)<(•))  M  (A  )  *((),$(• ))  i8  given  by 


(5.6) 


5  -  rA-Tx  and  ij>(s)  -  5  -  (s  +  Ox,  _r  <  x  i  °* 
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And,  from  Theorem  4.4,  P^(A^*)-*P^(0,4(O)  *  iNa^,  where 


a0  "  5 


(5.7)  a^k  =  5  +  l  (<jy)x  =  £  +  rx,  0  <  2k  <  N 

jodd 

l<j<2k-l 


a2k+l  5 


2N-1 


rx  +  2  l  (I)x 
jeven 
2£j<2k 


2N-2k-l 

C  -  (—  -4—  )rx,  0  <  2k  <  N. 

N  — 


Since  (*2k  +  a2k+l)/2  ~  *((tL  +  T2k+l)/2)  and  +  a^)/2 


iJ)((T2k_i  +  T2k^2^’  ^°^-^ows  from  (5.6)  and  (5.7)  that 


r 

"2JJ  2 

(5.8)  »PN(AN*r1PN(0,<^(.))  -  (A*)“l(0,*(O)«2-  2N  /  |2Nsx|2ds  -  L-  \x\2*  0. 

^  0  J 


Now  we  may  state  the  main  theorem. 


Theorea  5.4:  There  exists  an  element  z  €  Z  and  t  >  0  such  that 


SN*(t)z  does  not  converge  to  S*(t)z. 

Proof :  If  for  every  z  e  Z  and  t  >  0,  SN*(t)z  converges  to  S*(t)z, 

then  it  follows  from  Theorem  4.4  and  Lemma  5.2  that  for  every  z€Z 
PN(AN*)  *P^z  +  (A*)~*z,  where  by  Lemma  5.3  one  can  assume  that 
and  A  =  Aq  +  Aj  is  invertible.  This  contradicts  the  facts 
(5.8). 


A(.)  -  0 
(5.5)  and 
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